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(— I i Abstract. A microscopic description of Hawking radiation in sonic black holes has 

been recently presented [T| . This exactly solvable model is formulated in terms of one- 
dimensional scattering of a Fermi gas. In this paper, the model is extended to account 

C^ . possible finite size effects of a realistic geometry. The flovi' of particles is maintained by 

a piston (i.e. an impenetrable barrier) moving slowly towards the sonic horizon. Using 
, existing technologies the Hawking temperature can be of order of a few microkelvin in 

pH , a realistic experiment. 
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1. Introduction 

Analogue models for black holes |2j are recently subject of a renewed interest 
miaiSlilElinilZlElinillDllIIl since the realization of atomic Bose-Einstein 
condensate [121 CH! and Fermi degenerate gases [HI Unj at temperature of a few tens 
of nanokelvin. These achievements together with the realization of one- dimensional 
systems |16| [T7| ITHl ITUl (201 1^ make the experimental realization in the laboratory of an 
analogue of Hawking radiation [22] possible. This possibility has been made clear since 
the understanding of the relationship of sonic Hawking radiation with "the tunnelling 
through the top of the hill" 0. 

Analogue models for black holes have been originally considered in the hope that 
they can shed some light into the problem of quantization of gravity and on the related 
problem of information loss [2ni- In 1974, Hawking [21j proved that black holes can 
emit particles with a temperature proportional to its surface gravity. The predicted 
spectrum of the radiation is purely thermal and incoherent and this implies a loss of 
information after the black hole has evaporated away and disappeared completely. In 
1981, Unruh [221 suggested that an analogue of Hawking radiation can be produced in 
a quantum fluid turning supersonic. Very recently, a microscopic theory for this sonic 
analogue of Hawking radiation has been presented Ij. This theory, although unable to 
describe the process of evaporation of the black hole analogue, is based on elementary 
quantum mechanics and thus, it is unitary and information preserving. 

In this paper, I review the microscopic model for sonic Hawking radiation (Section 
2) and extend the model to a finite size geometry where the source of the particles is 
replaced by a moving piston, i.e. an impenetrable barrier (Section 3). This problem 
bears analogies with the "Fermi accelerator problem" [2^]- The possibility of an 
experimental realization of the model is discussed in Section 4. Session 5 concludes. 

2. Microscopic description 

A one-dimensional Fermi-degenerate non-interacting gas that scatters against a very 
smooth potential barrier provides a clear and simple quantum mechanical microscopic 
description of the sonic analogue of Hawking radiation. Sonic Hawking radiation is 
formed due to quantum tunnelling through the top of the potential barrier. Excitations, 
in a semiclassical sense are created at the barrier and are then radiated with a momentum 
distribution characterized by a temperature in exact agreement with Hawking-Unruh's 
formula fI^ |22l IHl- Although the evolution of the fluid's flow is unitary, when the 
momentum distribution is measured locally, for instance through the measurement of 
the dynamic structure factor, it is not distinguishable from a thermal (and incoherent) 
distribution. 



The sonic analogue of black hole radiation 3 

2.1. One- dimensional systems 

Consider a Fermi-degenerate gas of particles moving in a sufficiently narrow and 
smooth constriction such that their motion can be considered one-dimensional. Their 
stationary single-particle wave-functions can be written in the adiabatic approximation 
[SE] as '^{x,y,z) = ilj{x)^Q-x{y, z) where ^o.x{y,z) is the transverse ground-state wave- 
function with eigenenergy eo(a;). ^(x) satisfies a one-dimensional Schrodinger equation 
[—{h'^/2m)d'^ + Vex_t{x)]ip{x) = eip{x) where \4xt is an effective potential that includes 
also the zero point energy eo(x), i.e. 

Kxtlx) = (eO;x| - (^V2m)(5^2 + dl) + Kxt(x,l/,^)|eO;.) (l) 

where Kxt(a^, y, -z) is the three dimensional external potential. I assume Vcy±{x) very 
smooth and essentially non-zero only near its maximum Knax located at x = 0. 

2.2. Unitary description 

The transonic flow has an unitary quantum mechanical description. Consider a flow of 
non- interacting particles coming from the left and colliding against the potential Kxt- 
Their particle's wave-functions are given by scattering wave functions asymptotically 
defined by ip^ oc [e*'^^ -|-r(e)e~*'^^]6(— x) +t(e)e*''^9(x), where e = h^k"^ /2m is the energy 
eigenvalue and r(e) and t(e) are the reflection and transmission amplitudes, respectively. 
The evolution of the particle flow is unitary and the many-body wave function of the 
flow of particles is represented in second quantization by 

i/)= n «'!-) (2) 

0<e<emax 

where a\ is the creation operator of ip^^ and |— ) is the vaccum of particles. In this state 
all 'j/'e's are occupied up to an energy emax, which is assumed higher than V^ax- The 
non-interacting case is an ideal example of a fluid turning supersonic (see below) where 
the microscopic wave function of the flow (j21) is given from the beginning. 

2.3. Hydrodynamic description 

Since the barrier potential is smooth the particle propagation is semi classical and 
in ffist approximation |t(e)p ~ 0(e — Vmax)- Therefore the fluid maintains in a 
ffist approximation the sharp Fermi-degenerate character of the local momentum (or 
velocity) distribution, which represents the "Vacuum" of the black hole sonic analogue. 
This allows us to describe the fluid dynamic with the standard hydrodynamic equations. 
The velocities of the particles in the fluid are uniformly distributed within an interval 
{vl,vr). Figure 1 shows an example of velocity distribution. The right Fermi velocity 

vb^ = V^max - 2Kxt(a;)/m (3) 

is the velocity of a particle moving from the left to the right with initial velocity higher 
than the classical threshold velocity to go over the barrier 

Wesc = A/2Knax/"^ (4) 
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The left Fermi velocity vl is the velocity of a particle starting from the top of the barrier 
(located at a; = 0), i.e. 

Vl = sgn(x) ^/v^^^ - 2Vc^tix)/m. (5) 

The sgn(a;) function accounts for the two possible different signs of the velocity that the 
particle can have falling either to the right or left of the barrier. 




Figure 1. Semi classical velocity distribution as function of the coordinate x. 
Velocities are uniformly distributed in the area limited by vr (top solid line) and 
Vl (bottom solid line). Thinner solid lines correspond to different classical trajectories 
of the particles of the fluid. Here Emax = 2.5 and Vi^xt — exp(— x^) in dimensionless 
units. The inset shows the Mach number v/c as function of position x. 



The flow velocity v = {vl+vfi)/2 is the velocity of a local reference frame from where 
the fluid appears in equilibrium. Introducing the Fermi velocity vp = (— fL+'yR)/2 as the 
Fermi velocity in the co-moving frame, the density is given by n = rnvF/iTh. With these 
definitions the current Q = nv = rnvvF/irh is conserved. [The current conservation 
can be written in a form that is reminescent of the quantization of the conductance 
in one-dimensional channels: Q = {^r — fJ^L)/hp where hp is the Planck constant, 



Hl = "^Wesc/2 and fiR 



mvr, 



J'l\ The Bernoulli equation v^ jl + /i + Kxt/""^ = const 



where h[n) = 7r^/i^n^/2?T2^ is the enthalpy (the Fermi energy divided by the particle 
mass) is satisfied with const = f esc/2 -|- w^ax/2- The sound velocity (c^ = ndh{n)/dn) 
is consistently given by the Fermi velocity c = vp- It is easy to verify that the above 
defined particles flow is a transonic flow, namely subsonic on the left of the barrier and 
supersonic on the right. The inset of Figure 1 shows the behaviour of the Mach number 
v/c in the one-dimensional channel. 

Although, here I showed that the stationary hydrodynamic equations are satisfied 
for the gas fiow, it is important to state that also the long-wave perturbations of 
the stationary solution, i.e. sound waves satisfy the time- dependent hydrodynamic 
equations in the limit of very small amplitudes. This will be reported in details 
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elsewhere. This condition is important in order to apply Unruh's theory for sonic black 
holes 



2.4- Quantum effects: the temperature of sonic analogue of Hawking radiation 

Even if the potential barrier is very smooth, its finite thickness near its top smears 
the semi classical Fermi distributions at vi^{x). This is due to those particles colliding 
against the barrier with energies comparable with Knax- The reflection probability is 
given by (see also |i27^ ) 

\r(e)\^ = ^ f6) 

where Ux is the frequency of the inverted harmonic oscillator obtained expanding the 
barrier potential \4xt = Vmax — ^"m-^lx"^ + o{x^) around its maximum. Hawking's formula 
for a sonic black hole is given by |22l 0] 

^^ = (2^) -^ ' ^^^ 

where the surface gravity k in Hawking temperature Kh/2nkB is replaced by d{v — c)/dx. 
For the transonic flow of fermions defined above d{v — c)/dx = d{vL)/dx = ui^ and 
Hawking temperature ((Tj) becomes 

Therefore the particles reflected from the barrier are distributed according to the 
reflection probability ©, which is identical to a thermal Fermi distribution 

l + exp((e-/i)/fcBT) 
with a chemical potential fi = Knax and a temperature in exact agreement with Hawking 
temperature (jH)). 

When quantum tunneling is negligible the quantum state observed over a large 
segment L far away from the barrier would appear like a zero temperature Fermi- 
degenerate gas, which can be viewed as the semiclassical vacuum of an effective 
low-energy description. The role of the tunneling through the top of the barrier is 
to introduce a smoothed Fermi left point in L and thus create excitations on the 
semiclassical vacuum. 

We note that this interpretation of sonic Hawking radiation as related to a tunnel 
effect is consistent with the recent derivation of Hawking radiation as a tunnelling 
process of particles in a dynamical geometry due to Parikh and Wilczek [22] • 

Moreover as quantum tunnelling is very sensitive to the details of the potential 
barrier, the momentum distribution, related to the reflection coefficient of the barrier, 
can be different from the inverted parabola formula ©• This suggest that some 
information on the internal geometry of the black hole is carried outside by the radiation. 
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2.5. Measuring the sonic Hawking radiation 

Sonic Hawking radiation can be detected through the measurement of the dynamic 
structure factor S{q, u) 28| in an extended segment of the subsonic or supersonic region 
[T]. The S{q, u) of the flow is given P in the subsonic region by 

S{q,^) = SZ,{q,uj) e (^^) + 5^1, Jg,..) e (^-^) (10) 

and in the supersonic region by 

^(?,^)= K.(?,^) + ^tL,(?>-^)]0(^) (11) 

with /ii = Vjaax and fiR = Cmax- S^,ij.{1j^) is the dynamic structure factor of an 
equilibrium one-dimensional Fermi gas at temperature T and chemical potential fi given 
by 5^^(g, ^) = 1st e^/'-^/{Cosh{nnj/kBT) + Cosh[(^ + g - /^)/A:bT]} where L,s 
is an effective length of the probe-system interaction region [1 . Its T = limit is given 

by SZi<l, -) = tSi (-- + S + ^) (- + S - #) with A:^ = V2^,/h. 
Figure 2 shows some examples of S{q,Lj). S{q,uj) has a typical peaks structure at 
low momentum, which contrasts the non-interacting two and three dimensional Fermi 
gas 28 . The peaks are related to the dispersion relation u{q) = — wl \q\ Q{—q) + 
v-Ei\q\Q{q) valid for small q. The characteristic wave- vectors of Hawking radiation is 



qc = kBTu/hc ~ l/27r^/^?T, where Ix = y^m/hwx and n is the density. Indeed at very 
low frequencies huj <^ /cb^h and momentum q <^ qc the S{q, u) is dominated by the 
Hawking radiation. The peak values at \uj\ = \vi^q\ tends to 1/4 of the corresponding 
zero-temperature value L^sm / 2nh^q but with a half-width Au = 1.76 qkBTn/mvi^, 
which is much larger than the corresponding half-width Auj = hq^ /2m of the zero- 
temperature contribution. The measure of the width of these peaks provides a measure 
of the Hawking temperature of the sonic black hole. 

This calculation shows also that a local measurement of S{q, uj) through a scattering 
experiment in the subsonic or supersonic regions does not probe the underline global 
coherence of the flow. The S{q,uj) in the subsonic region is equal to that of a clean 
one-dimensional channel connecting two reservoirs: one on the left at zero temperature 
and one on the right at the Hawking temperature Th given by (jH))- This represents a 
one-dimensional analogue of Hawking radiation from a black hole into the surrounding 
zero-temperature space. 

To distinguish between the pure quantum state of this fluid, which is represented by 
the many-particle wave function Q and a mixed thermal state using only measurements 
made in the subsonic or in the supersonic region is probable difficult. This would require 
the measurement of correlations such as those between opposite sides of the barrier. 
When a particle is detected in the subsonic region at some negative momentum /c, i.e. 
reflected by the barrier, we know that it will not be detected in the supersonic region 
at momentum —k. 
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Figure 2. Sonic Hawking radiation can be detected through the measurement of 
S{q,uj). The figure shows the w-dependence of S{q,uj) in the subsonic (left) and 
supersonic (right) regions and for a positive (top) and a negative (bottom) value of 
q. Here fceTH — 0.15 Vmax, Cmax — 2.5 Vmax and \q\ = 0.05mVesc/h is chosen in the 
region of small q where thermal-like transitions (broadened peaks) dominate over the 
zero-temperature contributions (step-like peaks). S{q,uj) is in units of Lcff'm/2Trh'^q 
and Lo is in units of Knax/^- 



3. The piston geometry 

In the model introduced in Section 2 the quantum mechanical description of the many- 
body system is given in terms of single-particle scattering wave-functions. This is a very 
ideal situation where we theoretically identified the predicted Hawking radiation. In a 
practical realization of the model, one can perhaps imagine to join the one- dimensional 
channel with a large reservoir on one side the channel to ensure a continuum source for 
degenerate Fermi particles and with an empty area on the other side where the atoms 
can escape. One needs to assume that in the reservoir the interaction of the particles 
is sufficiently effective to ensure that the quantum mechanical statistical equilibrium 
takes place quickly enough. Therefore, the complete quantum mechanical description 
of the entire setup is not possible in simple terms and the local quantum mechanical 
description in the one-dimensional channel relies on assumptions on the reservoir. 

In the following we consider a finite size system as a model for a sonic black hole 
that can be fully described quantum mechanically and no assumption of equilibration is 
necessary. The reservoir is replaced by a time- dependent impenetrable potential barrier 
(the "piston" ) that with his constant motion pushes the gas particles from the subsonic 
region over a smooth barrier in the supersonic region. The transonic flow will take place 
at the top of the smooth barrier. 
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3.1. Hydrodynamic description 

The piston, i.e. an impenetrable and sharp potential barrier of coordinate Xpistonit) < 
moves at constant velocity fpiston > in a region of zero potential towards the barrier 
located at x = 0. The Fermi-degenerate gas is initially trapped between the piston and 
barrier, i.e. in the subsonic region. The piston compresses the gas, creating thus the 
conditions of pressure to start the transonic flow. The gas will later flow stationary for 
a long time and its velocity v in the subsonic region will be equal to the piston velocity, 

i.e. V = fpiston- 

In the initial part of the gas dynamic, i.e. before that the atoms start to escape to 
the right of the barrier, the gas particles have an uniform semiclassical distribution of 

velocity in the subsonic region given approximatively by 

m 
f{x, V, t) dxdv = —- Q[—vi,{x, t) — v] 0[t>pj,(x, t) — v] dxdv (12) 

Alp 

where h-p is the Planck constant. Initially the local left and right Fermi velocities vi,{x, t) 
and fR(a:, t) have modulus less than the escape velocity Wesc and are given by 

X 

I •''piston l^'^J I 
X 

|'^pistonl,''j I 

where 



t;L(a;,t) = -t;F(t) -Wpiston-r: 7ITT (^3) 

V^{X, t) = VF{t) - ^'pistonT— ^-TITT (14) 



•^piston l^f-j -^piston l,Uj + fpiston f- \^^) 

The Fermi velocity VF{t) is initially fixed by the conservation of the number of fermions 
Np in the left side region and is given approximatively by 

m\Xpiston{t)\ 

We assume for simplicity that the barrier size is much smaller than |a;piston(^)|- The finite 
size of the barrier results in a modification of the transient dynamics, which discussion 
is beyond the purpose of this paper, and does not concern the quantum effect occurring 
at the sonic event horizon. 

Every time that one atom bounces against the moving piston its absolute value of 
the velocity increases by 2t'piston- At a time ttrans some of the particles near the barrier 
reach the escape velocity fesc and they start to go over the barrier in the supersonic 
region. After a transient time of about 2 1 Xpiston (i^trans ) |/f^esc the gas will enter a stationary 
regime characterized by a semiclassical distribution of velocity with a form still given 
by equation (fT^ but with constant left and right Fermi velocities vp and wr given in 
the subsonic region by 

f L = - ^^esc (17) 

Vr = t^esc + 2t;piston (18) 

and in the supersonic region by 

Vl = ^^esc (19) 

^^R = ^^esc + 2t;piston (20) 
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Figure 3 sketches the seiniclassical Fermi distribution during the stationary phase of the 
gas dynamics. Note that the particles that colhde with the piston at velocity wl are 
reflected back at velocity fR. Thus the moving piston connects v-^ with v^.- 
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Figure 3. Semi-classical velocity distribution as function of the coordinate x for a 
certain time chosen in the stationary regime. Velocities are uniformly distributed in 
the area limited by the solid line. Thinner solid lines correspond to different classical 
trajectories of the particles of the fluid. Here T^xt — (1/2) cxp(— a;^), Ucsc = 1 and 
Vpiston = 0.25 in dimcnsionlcss units. Note that the piston (x « —12) moving to the 
right recombincs the left Fermi velocity v^ (bottom part of the solid line) with the right 
Fermi velocity vj^ (top part of the solid line) in the subsonic region. [For the purpose 
of rendering the particles trajectories in this "instantaneous" picture it is necessary to 
assume that the collisions of the particles with the piston take place instantaneously.] 



The gas flow will eventually terminate when the piston reaches the barrier after 
the time interval |a;piston(^trans)|/i'piston- The above stationary transonic flow can also 
form starting with different initial conditions and these will affect mainly the transient 
dynamics. For instance, the transonic regime can start also by lowering the barrier high. 



3.2. Quantum effects of the sonic event horizon 

The next step is to consider the quantum description of the gas dynamics and the 
consequences of the tunneling at the top of the barrier. The quantum state of the Np 
particles gas is given by 



I/) 



n 



-) 



(21) 



l<n<Np 



where |— ) is the vaccum of particles and a|^ creates one fermion in the single-particle 
state ipnixji), which is solution of the time- dependent Schrodinger equation. In the 
initial stage of the dynamics ipai^jt) are quasi eigenstates (resonant state) of the 
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subsonic region multiplied by a common phase factor exp [iS{x, t)\] that fixes the initial 
hydrodynamic velocity field. In the limit of a negligible barrier size the time- dependent 
potential is equivalent to the boundary conditions ^'(^^piston,^) = i'{0,t) = 0. Therefore 
the ip^{x,t) can be approximated in the initial stage of the dynamics by 



tpn{x,t) = sin 

V 2 ppistonl^jl 



nil X 



I '^piston l^t-j I _ 



exp 



2 



— 1 



''^ '^piston X 
Z,n \Xpiston\J^ ) I 



(22) 



which are exact solutions of the Schrodinger equation with the above specified boundary 
conditions. Note that as the piston evolves towards the barrier the number of nodes n 
of V'n is conserved in time and its effective wavelength grows. The semiclassical velocity 
distribution of (PT|) with the wave-functions ipn{x,t) given by (j^ corresponds indeed 
to the initial distribution given previously by equation (jT^ and equations (fT^. (fT^ . 
(fTHjl and (fTHjl as one can easily check. 

As the piston moves even closer to the barrier the wave- functions ipnixjt) start to 
propagate one by one over the barrier. The differences in energy levels in the subsonic 
region are assumed to be much less that the Hawking energy, i.e. /ifesc/2|xpiston| *^ ^b ^h 
where Th is given by equation (jH)). We can thus use the transmission formula (jHl) to 
predict the distribution of momenta in the supersonic region. The predictions will be 
compared with the results obtained by the numerical integration of the Schrodinger 
equation in subsection 3.3. 

Consider a particle transmitted in the supersonic region with some velocity f final- 
Before exit the subsonic region it has been reflected between the piston and the barrier 
a certain number of times. Therefore the probability p of this process to occur is 
proportional to 

P(^final) OC |t(t'final)r ^(^final " 2i;piston)P ••• ^(^final " 2 TlcoU t^piston) ^ 

(X |t(t;finai)|' n:^='i'|r(t;finai - 2nt;pi,ton)|' (23) 

where ricoii is the number of collisions with the barrier and |tp and |rp are the 
transmission and reflection probabilities as function of the velocity. In the stationary 
regime, the flnal distribution of velocity f{v) should not depend anymore on the initial 
distribution and it should be given, a part from a normalization constant, by 

fiv) = \tiv)\'lC=T\riv - 2nt;pi3ton)r (24) 

where nupper is an integer for which |r(t;esc - 2raupper fpiston)^ = 1- 

In particular, when the Hawking energy is much smaller than the difference between 
left and right Fermi energies in the stationary regime, i.e. /cb ^h <^ ^ "ypiston "^esc equation 
()24|) simplifles to 

fiv) = \tiv)\^\r{v-2vpi,ton)\^ (25) 

as the next terms of the product in P^ are all equal to one. This condition is equivalent 
to the condition k^ Th -C mc'^ where c is the speed of sound on the event horizon. Thus 
in the supersonic region, the occupation probability p(e) of an energy state e is given a 
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part from a normalization factor, by 

„/^N exp((e - K,ax)/fcBT^) ,261 

^^' [1 + exp((e - K,ax)/fcBTO)] [1 + exp((e' - V^,^)/kBT^)] ^ ' 

where e' is the energy of the particle before colliding to the piston 

e'(e) = m( v/20^ - 2t;piston)V2 (27) 

Therefore differently from the case discussed in Section 2 both the left and right Fermi 
points vi,{x) and fR(x) are smoothened by the tunnelling through the barrier. The 
reason of this is that wl and wr are recombined by the piston (see Figure 3). Remarkable 
the gas appears as thermally distributed on both sides of the local Fermi velocity 
distribution. 

3.3. Numerical results 

We would like to compare the above theory with a numerical integration of the 
Schrodinger equation for the many-particle wave- function of the gas (j21|) . However, 
even if the system is non-interacting and therefore treatable in principle we should note 
that the above conclusions are valid for a very large system. Different requirements need 
to be simultaneously satisfied: (i) smooth and relatively extended barrier; (ii) IxpistonI 
much larger than the barrier extension to minimize finite size and transient effects; (iii) 
large number of fermions; (iv) long time scale as the piston moves slowly towards the 
barrier. 

On the other hand, the time evolutions of the wave-functions ipn are quite similar. 
Each -011 defined in (j^ is a superposition of exp(i nx / |xpiston| —irn f piston x'^/2h |xpiston|) 
and exp(— iraa;/ |xpiston| — i^^'^piston a;^/2/i IxpistonI), which contain classical velocities 
uniformly distributed in the intervals (/in/m |xpiston|, ^'^/'"^ |a^piston| + "ypiston) and 
(— /in/m IxpistonI, — ^^/'"^ l^^pistonl + "^piston), respectively, where h is the Planck 
constant. Each ip^ performs thus a reasonable sampling of the phase space distribution 
of the gas. We therefore study the time-evolution of only one of the ip^- 

We calculate the stationary momentum distribution of the flow ()21|) . which is a 
Slater determinant by approximating the sum of single-particle contributions by only 
the contribution of one of the ipn- Its momentum distribution in the supersonic region 
is calculated from the frequency components of ^n(^ref5^) at some arbitrary reference 
point Ti-ef in the supersonic region where the particles propagate as a plane waves and 
taking into account the density of states. 

Figure 4 shows some examples of the comparison between the momentum 
distribution in the supersonic region calculated from one of the ip^ and equation ()24|) with 
(jHl) as formula for the reflection probability. The agreement is reasonable considering 
that the chosen initial wave functions have only n = 14 nodes, ip^ is a quasi-resonant 
state of the left region multiplied by the phase factor exp [— imt;piston x^/2h |xpiston(0)|], 
which gives a constant gradient of velocity field decreasing from the maximum value 
I'i^pistonl at the classical turning point of ip^^ at Xpiston to zero near the barrier. 
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Figure 4. (a) Example of numerically computed density of an initial resonant state 
with n = lA nodes (solid line) and potential energy (dashed line) at the initial time of 
the evolution. In dimensionless units Knax = 1/2, Vesc = 1 and fcsT]^ — 0.02. Figures 
(b), (c), (d) show the occupation probabilities in the supersonic region of plane waves 
with energy e obtained by numerical integration of the Schrodinger equation (solid 
lines) and from the analytical formulas (|24|l and ^ (dashed lines). In figure (b) and 
(d) (fcsT^ = 0.02 and fcsT^ = 0.04, respectively) Upiston = 0.18 and the dashed hne 
is calculated using 1)25(1 where the multiple reflections corrections (|24|l are neglected. 
In figure (c) Wpiston — 0.015 is chosen such that the Hawking energy ksT^ = 0.02 is 
comparable to mc^ at the horizon and the multiple reflections corrections 124|) (dashed 
line) are necessary. Dot-dashed line is obtained from (|25|1 . Note that in all three 
examples (b), (c), (d) the classical gas distribution is always a sharply defined interval 
showed for comparison with the dot-dashed line only in (c). The small deviations 
between solid and dashed lines are quite sensitive on small changes in the exponential 
phase factor (not shown here) and on the number of nodes of the initial wave function, 
which is n = 14 in all three examples. 



4. Hawking radiation in the laboratory: Relevant experimental parameters 



Today's atomic trapping techniques are sufficiently developed to ensure the realization 
of a setup for the demonstration in the laboratory of the sonic analogue of Hawking 
radiation. The main difficulty is the realization of an "event horizon" barrier that allows 
appreciable tunnelling. Josephson-like effects in Bose-Einstein condensates [301 1^ have 
been recently observed jH2]- This proves that it is possible to achieve barriers that 
allow a reasonable tunnelling. Here follows an estimation for a Fermi degenerate gas of 
Lithium atoms, which has been recently obtained jBH!- Lithium can have a degeneracy 
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temperature as high as 8 /iK, which fixes the scale of emax- The barrier potential could be 
created by a far-blue-detuned laser sheet propagating in a direction orthogonal to that 
of the fiow propagation. One of the resonance lines of Lithium is at 233 nm [HI] . Using 
an ultra-violet diode laser, which is commercially available and tightly focused with a 
waist (1/e^) of twice the laser wavelength (that can be tuned at 233 nm) the Hawking 
temperature is of order of Th ~ 2 /xK, which is a reasonable value and comparable 
with the degeneracy temperature of 8 /xK. The characteristic wavelength of Hawking 
radiation Ac is of order of 4 /xm. 

5. Conclusions 

The exactly solvable model, originally formulated in terms of one-dimensional scattering 
of a Fermi gas p] is extended to account for a finite size realistic geometry. The fiow 
of particles is maintained by a piston moving slowly towards the sonic horizon. Both 
points of the Fermi segment appear thermally smoothed. Moreover multiple refiections 
between piston and barrier lead to possible modifications of the result obtained in P] 
for the momentum distribution. Interestingly, due to the geometry of the moving piston 
the prediction on the thermal momentum distribution can be tested studying the time- 
evolution of only one single-particle wave-function. 

We estimated that by using existing technologies the Hawking temperature can be 
of order of a few microkelvin in a realistic experiment, which is a reasonable temperature 
to be measured. 
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